Berry phase, semiclassical quantization and Landau levels 
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We propose the semiclassical quantization for complicated electron systems governed by a many- 
band Hamiltonian. An explicit analytical expression of the corresponding Berry phase is derived. 
This impact allows us to evaluate the Landau magnetic levels when the rigorous quantization fails, 
for instance, for bilayer graphene and graphite with the trigonal warping. We find that the magnetic 
breakdown can be observed for the certain type of classical electron orbits. 
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The most accurate investigation of the band structure 
of metals and semiconductors is studying the Landau 
levels in magneto-transport and magneto-optical exper- 
iments. However, the theoretical solution of the band 
problem in magnetic fields cannot often be exactly found. 
A typical example is presented by graphene layers. For 
bilayer graphene and graphite, the effective Hamiltonian 
is a 4 X 4 matrix giving four energy bands. Fig. [T] shows 
nearest two bands of the level structure together with 
semiclassical orbits. The trigonal warping described by 
the effective Hamiltonian with a relatively small param- 
eter 73 provides an evident effect (see right panel). An- 
other important parameter is the gate-tunable bandgap 
U in bilayer graphene. In this situation, the quantization 
problem cannot be solved within a rigorous method. To 
overcome this difficulty one can use a perturbation the- 
ory, however this theory becomes quite complicated for 
the many-band Hamiltonian. 

Alternatively, the semiclassical quantization can be ap- 
plied. Thus, we can use the Bohr-Zommerfeld condition 
as 
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Here iS'(e) is the cross-section area of the electron orbit in 
the k space for the energy e in absence of the magnetic 
field B and for the constant momentum projection k^ 
on the magnetic field, n is an integer supposed to be 
large. T is the number of the smooth turning points on 
the electron orbit. There are two smooth turning points 
for the Landau levels and only one for skipping electrons 
reflected by the hard edge. 

The goal of this letter is an explicit analytical expres- 
sion for the (5(e)— phase within the band scheme of the 
matrix Hamiltonian. The semiclassical approach is used 
for the magnetic field normal to the layered system when 
the quantization of in-layer momentum components is 
only essential and the size of the Fermi surface is small 
compared with the Brillouin zone size. We illustrate our 
results for bilayer graphene. Notice, that the (5(e)— phase 
depends on the energy and can be taken in the inter- 
val < |(5| < 1/2. If the spin is neglected, (5 = and 
T = 2 for the Landau levels, and 5 = 1/2 and 7" = 2 



for monolayer graphene. In these two cases, the semi- 
classical result coincides with the rigorous quantization 
and it is closely connected with the topological Berry 
phase ll|. This (5— phase was evaluated for bismuth in 
Ref. [21, preceding Berry's work by almost two decades, 
and it was considered again for bismuth in Ref. 13|. For 
graphite, the semiclassical quantization was applied in 
Ref. [J|. However, in the general case, the evaluation 
of the (5— phase is still attracted a widespread interest 

M- 

The problem under consideration is described by the 

Hamiltonian in the band representation 



(V-k + r-£)^ = 0, 



(2) 



where the column '5 consists of functions corresponding 
with a number of bands included and is labelled by the 
band subscript which we omit together with the matrix 
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FIG. 1: (Color online) (a) The energy dispersion e{k, a) of two 
nearest bands (the electron band shown in solid line and the 
hole band in dashed line) in bilayer graphene for two polar 
angles a with the local extrema at fc 7^ ("mexican hat") 
represented. The band parameters are given in the figure, 
others are 70 = 3.05 eV, 71 = 360 meV, 74 = —150 meV 
[12. I13I]. (b) Cross-sections k{a,e) of the electron band for 
energies of 80 meV (dashed-dotted line) and 40 meV (solid 
line) . 



subscripts on T and V; a summation over them is implied 
in Eq. ([2]). Matrices F and V are the first two terms in 
a series expansion of the Hamiltonian in the power of 
quasi-momentum k. 

In the magnetic field, the momentum operator k de- 
pends on the vector-potential A by means the Peierls 
substitution, 

k — —ifiV - eA/c, 

providing the gauge invariance of the theory. The mag- 
netic field can also enter explicitly describing the mag- 
netic interaction with the spin of a particle. However, 
for the graphene family, the magnetic interaction is weak 
and omitted here. 

A simple example of Eq. ([2]) is given by the graphene 
monolayer. There are two sublattices in it, and Eq. ([2]) 
is represented by a 2 x 2 matrix if the spin of carriers is 
neglected. Another example considered below is bilayer 
graphene with the 4x4 matrix Hamiltonian. For the 
monolayer and bilayer graphene, both V and k are two- 
dimensional vectors, e.g., with x and y components. 

We seek for \1/ in the form 



vj/ := $ exp {is/h) , 



(3) 



where the function s is assumed to be common for all 
components of the column \l/. The equation for $ is 
reduced to 

[V-(k-inV)+F-£]$ = 0, (4) 

with k = Vs - eA/c . (5) 

The function $ is expanded in series of Ti/i: 



* = V f - I 'Pn 



Comparing the terms involving the same powers of Ti in 
Eq. (HI we have 



(V • k + F - e)^™ = -VV(^„,_i 



(6) 



For TO = 0, we get a homogeneous system of algebraic 
equations 



(V-k + F-e)^o = 
which has a solution under the condition 
Det(V-k + F-e) = 0. 



(7) 



(8) 



This equation determines the classical electron orbit, 
e{kx, ky) = e, in presence of the magnetic field while the 
electron energy s is constant. At the same time, the 
equation yields the electron dispersion equation with k 
as the momentum without any magnetic field. In 3d case, 
the electron dispersion depends as well on the momentum 



projection kz on the magnetic field and our scheme can 
be implied in this case without the expansion in kz- 

It is convenient to choose the vector-potential in the 
Landau gauge A^ = —By, Ay = Az — in such 
a way that the Hamiltonian does not depend on the 
x— coordinate. Then, the a;— momentum component K^ 
becomes a conserved quantum number and the function 
s in Eq. ^ can be written as 



s = xK,j; + a{y) . 
The equations ^ are reduced to 



(9) 



kx = Kx + -By , ky ^ -—. 
c dy 

These equations enable us to use the variable k^ instead 
of y and to obtain 



c 
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<^(^^) = 7^ / hiK)dk'x 



(10) 



where ky as a function of kx is determined by the disper- 
sion equation dH). 

The eigenfunction column obeying Eq. ([7|) can be mul- 
tiplied by the scalar function C common for all elements 
of the column 

ifo -^ Cipo 

where ipo is any eigen-column of Eq. ([7]). The function C 
is determined by Eq. ^ with m = I. Left-to- right mul- 
tiplying both sides of this equation by ipQ and using the 
Hamiltonian hermiticity, i.e. the complex conjugations 
of Eq. ([7|), we get the consistency condition 



^SV.V(C^o)=0, 



(11) 



where the derivative with respect to y (i.e. to kx) is only 
to be taken. The left hand-side of this equation can be 
written as 



1 dC 
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C dkx ^ 2ip*Vyipa dkx ^ iplVyipo ° ^ dkx 

Using the identity (p^VipQ — (fiQipov with the electron 
velocity v — de/dVi, one can write the solution of Eq. 
dnj as 



C = co{ipa(poVy) ^/^ exp(-i6 



where 



= Im 



d(po 



dkx ,,, 

~1^ 'Po^y-JT 

iPofoVy dk. 
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and Co is the normalization factor. 

The quantization condition can be written as usual 
from the requirement that the wave function has to be 
single- valued. Continuing Eqs. (ITUl) . P^. and P^ along 



the orbit and making the bypass in the complex plane 
around the turning points where Vy = to obtain the 
decreasing solutions in the classically unaccessible region, 
one obtains 7~ = 2 and (5— phase as a contour integral 
along the classical orbit 
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Using the Hamiltonian hermiticity, after the simple alge- 
bra (see Ref. [4I), Eq. P^ can be rewrite as 



(5(e) = — Im (f> 

^ ' An J ip*ipo 
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called usually the Berry phase. 

Now let us calculate the 5— phase for bilayer graphene. 
In simplest case, the effective Hamiltonian can be written 
(see, for instance Refs. [l2,ll3l) ^^ 
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where the parameter U describes the tunable gap, 71 is 
the nearest-neighbor hopping integral energy, the matrix 
elements are expanded in the momentum k± = ^ik^ — ky 
near the K points of the Brillouin zone, and the constant 
velocity parameter v is incorporated in the notation q± ~ 
vk±. 

Here, the orbit is the circle defined by Eq. ([8|), written 
in the following form 

[([/ + ef - q\{U - ef - q'] - 7?(£2 - C/2) - . (16) 

The eigenfunction ipo of the Hamiltonian (jlSp can be 
taken as 



<Po 



f{U-s)[{e + Ur-q']\ 
q^q^-is + Uf] 

liq+{U ~e) 



V 



(17) 



with the norm squared 



vl^o = [{e + Uf-q^f[{e~Uf + q^ 



(18) 



The derivatives for Eq. ([H|) are calculated along the 
trajectory where the energy e and consequently the tra- 
jectory radius q are constant. The equation (J16l ) has 
only one solution for g^ if |t^l < kl < y/U"^ +7i- First, 
let us consider this case. 

(i) there is only one orbit at given energy e with the 
radius squared 



U' 



+ ./4C/2£2 



u'H 



The matrix Vy = dH/dky in Eq. p4)) has only four 
nonzero elements Vyi2 — Vy2i — Vyn = K/43 = — 1- 
Using Eqs. ^ and ^, we find 



dm 
' dk.^ 



Im vlVy^ = AUe{U - e)[{e -f Uf ~ q^] 



(19) 



This expression is constant on the trajectory as well as 
(pl^Po, Eq. dH]). Therefore, in order to find 5, Eq. HH), 
we have to integrate along the trajectory 



This integral equals —dS{e)/de, where ^(e) = nq'^ is the 
cross-section area, Eq. ([T]), with 



dS{e) _ 2(g- 
de 



C/2 - £2 



-7? 



g2-[/2 



(20) 



Now we have to substitute Eqs. UB^ UM^ and ([20 
into Eq. ^^. Thus, we find the Berry phase 



Sis) 



-eU 



-eU 



g2 - £2 _ JJ2 



V4C/2e2 + (e2 _ f/2)^2 



(21) 



shown in Fig. [2l where (5— phase of bilayer graphene with 
trigonal warping is also shown, the detailed calculations 
will be elsewhere published. For the ungaped bilayer, 
U = 0, the Berry phase ^(e) = 0. The Berry phase 
depends on the energy and S — =Fl/2 at e = ±U. At the 
larger energy, e ^ U, the Berry phase 6 — )■ =Ff^/7i- 

Substituting Eq. ()2ip in the semiclassical quantization 
condition, Eq. ([T|), and solving the equation obtained 
for e, we get energy levels as a function of the magnetic 
field. We have to notice that the Landau numbers ul 
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FIG. 2; (Color online) Semiclassical phase vs energy in the 
conduction band of bilayer graphene without trigonal warping 
(solid line) and with warping (dashed line). 



listed in Fig. |3] do not coincide with the numbers n in 
the semiclassical condition ([l}. The rigorous quantiza- 
tion shows that there are only one Landau level with 
riL = and three Landau levels with n^ = 1 [14]. These 
levels are not correctly described within the semiclassical 
approach. However, for ul > 2, there are levels in all 
four bands s (two nearest bands with s = 2, 3 are shown 
in Fig. [3]). They correspond with the quantum number 
n — ul — 1, and the semiclassical levels become in excel- 
lent agreement with the rigorous solution for the larger 
n. 

(ii) for |L/|/Vl + (2C//7i)2 < |e| < \U\ , at the given 
energy, there are two orbits with the radius squared 
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This is an effect of " the mexican hat" . Then we seek for 
the general solution as a sum of two solutions (ySgCl) ^nd 
(y5Q(2) corresponding to these two contours, 

with two scalars Ci and C2. Instead of Eq. ([Tl]) we have 
a system of two equations written in the 2x2 matrix 
form as follows 

,22, 



a 



dqx 



bC = 



where the notations of the matrix elements are intro- 
duced 



atk = ^oi'^)Vyipoik) , hik = ^oi'>-)^y 



dqx 
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FIG. 3: (Color online) Energy levels Esnj^ for the K valley 
in magnetic fields for bilayer graphene within rigorous quan- 
tization (solid lines) and in the semiclassical approximation 
(dashed-dotted lines); in the notation {sul), ul is the Lan- 
dau number and a — 1,2,3,4 is the band number, only two 
nearest bands (s — 2, 3) are shown at given ul from 2 to 7. 
There is only one level, |10), with ul ~ and three levels 
(s = 1,2,3) with TiL ~ 1. The levels for the K' valley are 
obtained by mirror reflection with respect to the £ = axis. 



The off-diagonal matrix elements aik vanish, ai2 = a2i — 
0. Thus, the first equation of the system (|22|) becomes 

dC 
2qiyr—-^ + (2*£^ - rqx/qiy)Ci + i{2eU + r)C2 = , 

^qx 

and the second equation can be obtained with the index 
replacement 1 o 2 and r ^>- —r . 

These equations can be simplified with the substitution 

C. = a(gf-g2)-i/4^ (23) 

For the new functions C.^, we get the equation system 



dCi 



l.^ 



dqx V ?2iy 2 

q2y- iEC2 -I. —^{E - -)Ci =0, 

dqx \l qiy 2 



where the parameter qiy = y^qf — q^, i = 1,2 and E = 
eU/r. 

For the minimum of conduction band (maximum of 
valence band), where r — >■ 0, there is a simple limit. 



dCi I ~ ^ . , ~ 

qiy- -61 = with 62 

dqx 2 

Solving this equation, one gets 



Ci — Co exp [ — arcsm — 
2 qi 



'Ci 



(24) 



Going with qx along the trajectories and making the 
bypass in the complex plane around the turning points 
qx — ±qi and qx — ±^2, we see that both Ci and C2 ac- 
quire from two turning points in Eq. (j23p the additional 
phase — TT with T = 2. At the same time, we have —1/2 
from Eq. (p4)) for (5— phase. Thus, at the boundaries of 
the narrow interval considered, the 5— phase obtains the 
same value, S = —1/2. 

Taking into account the phases of the functions ^p^ii), 
we see, that the area S{e) in Eq. ([T]) can play the differ- 
ent role. In weak magnetic fields, slower oscillations with 
the smaller S{e) corresponding to (72 should be observed 
in oscillating phenomena. However, when the magnetic 
field becomes larger and the semiclassical condition is 
fulfilled only for the larger cross-section ^(e), calculated 
with qi, the larger frequency oscillations should be ob- 
served. This is nothing but the magnetic breakdown ^15| 
which should be utilized if the chemical potential belongs 
to the interval where the effect of " the mexican hat" ap- 
pears. 

In conclusion, our study shows that the semiclassical 
approach gives a powerful tool for probing the electron 
magnetic properties in metals. The Berry phase depend- 
ing on the energy can be calculated and observed even for 
complicated band scheme. The method presented here 
should be useful for many electron systems. 
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